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ABSTRACT 

We present Sunyaev-Zel'dovich Effect (SZE) scaling relations for 38 massive 
galaxy clusters at redshifts 0.14 < z < 0.89, observed with both the Chandra 
X-ray Observatory and the centimeter-wave SZE imaging system at the BIMA 
and OVRO interferometric arrays. An isothermal /5-model with central 100 kpc 
excluded from the X-ray data is used to model the intracluster medium and 
to measure global cluster properties. For each cluster, we measure the X-ray 
spectroscopic temperature, SZE gas mass, total mass and integrated Compton-?/ 
parameters within r 2 5oo- Our measurements are in agreement with the expec- 
tations based on a simple self-similar model of cluster formation and evolution. 
We compare the cluster properties derived from our SZE observations with and 
without Chandra spatial and spectral information and find them to be in good 
agreement. We compare our results with cosmo logical numerical simulations, and 
find that simulations that include radiative cooling, star formation and feedback 
match well both the slope and normalization of our SZE scaling relations. 

Subject headings: galaxies: clusters 



1. Introduction 

The Sunyaev-Zer dovich Effect (SZE) is a unique and powerful observational tool for cos- 



mology (for review see lCarlstrom et al.ll2002l ). It is a small distortion in the cosmic microwave 
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background (CMB) spectrum caused by scattering of CMB phot ons off a distribution of high 



energ y electrons in dense structures such as clusters of galaxies (jSunyaev fc Zerdovichlll970 



19721 ) . This effect has a unique property that the signal i s independent of redshift, mak- 



ing it particularly well suited for deep cluster surveys (e.g., iHolder et al.lhoool: Iwel 



2004 



2002) . Several SZE survey experiments are currently in progress ( Ruhl et al.l 12004 ; iFowler 



or et al. 



Kanekoll2006l ). and are expected to generate a large sample of SZ-selected clusters with 
masses greater than ~ 2x 10 14 M Q . The resulting large samples of galaxy clusters will enable 
direct measurements of the evolution of the number density of galaxy clusters as a function 
of redshift and in principl e can provide a power f ul constraint on the nature of dark en ergy 
JWang fc Steinhardtl[l998l : Iviana fc Liddl3ll999l : iMohr et al.lboool : Irlaiman et allboOlh . 



To utilize the upcoming SZE cluster surveys for cosmological studies, it is important to 
understand the relation between the SZE observables and the mass of a cluster. If the evolu- 
tion of clusters is dominated by gravitatio nal processes, a simple model of cluster formation 
and evolution based on the virial theorem (lKaiserlll986l ) predicts simple power-law relations 
between cluster masses and certain integrated cluster properties, including the integrated 
SZE flux (which is proportional to Y, the integral of the Compton-y parameter over the 
solid angle of the cluster). Numerical simulations further suggest that Y sh ould be an ex 



2004; 


Motl et al. 


2005; 


Naeai 


2006) 



[e.g., 



da Silva et al. 



shift evolution of the SZE scaling relations are relatively insensitive to the details of cluster 
physics, although numerical simulation s show that the input cluster physics affects the nor- 



malization of the SZE scaling relations (lNagaill2006l ). It is therefore important to investigate 



the properties of the SZE scaling relations observationally. 

Previous studies have addressed the correlation between the SZE signal and X-ray prop- 
erties. For instance, ICoorayl (119991 ) found a positive correlation between the central SZE 



decrem ent and the X-ray luminosity in a sample of 14 clusters. Similarly, McCarthy et al. 



( 120031 ) detected correlations between the central SZE decr ement and X-ray dete rmined mass, 
temperature, and luminosity for a 22 cluster sample, and iMorandi et al.l (120071 ) for a sample 
of 24 clusters. These studies use data from multiple SZE and X-ray experiments, making 
systematics more difficult to control, and focus on the r elationship be t ween the central values 
of the SZE signal with the X-ray properties. Recently, iBenson et al.l (120041 ) showed that the 
integrated SZE flux is a more robust observable than the central values of the SZE signal, and 
found a strong correlation with X-ray temperatures usi ng a sample of 15 clusters obtained 



by th e Sunyaev-Zeldovich Imaging Experiment (SuZIE. iHolzapfel et al.l 119971 ; IBenson et al. 
20031 ) and X-ray temperatures from the ASCA experiment. 



This paper is the third in a series of papers combining SZE and Chandra X-ray measure- 
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ments of galaxy cluster s to study cosm o logica l properties, following iBonamente et al.l (120061 ) 
(B2006 hereafter) and lLaRoque et al.l (120061 ) (L2006 hereafter). Here we present observa- 
tional studies of SZE scaling relations for clusters of galaxies. This paper advances the results 
of previous cluster scaling relation works in several ways. First, we use the largest observa- 
tional sample yet constructed (38 clusters at redshift 2=0.14 — 0.89). Second, our analysis is 
based on SZE and X-ray observations obtained using the same instruments: all SZE fluxes are 
determined using centimeter- wave inter ferometric data from the BIMA/OVRO SZE imaging 
experiments (e.g. lLaRoque et al.ll2003l ). and all cluster X-ray properties are derived using 
data from the Chandra X-ray Observatory. Finally, our Chandra observations have an order 
of magnitude better spatial resolution than the X-ray data used in previous studies, which 
greatly improves our ability to identify and exclude compact foreground sources which are 
superimposed on the cluster X-ray emission. 

Throughout the paper, we assume a ACDM cosmology with Qm=0.3, ^a=0.7 and 
h = 0.7, where h is defined such that Hq = 100 h km s _1 Mpc -1 . All uncertainties are at 
the 68.3% confidence level. 



2. Theory of cluster scaling relations 

2.1. The virial radius and r 2 5oo 

In order to establish relationships between mass, SZE flux and other cluster properties, 
one needs to define a radius out to which all quantities will be calculated. This radius 
should be physically motivated, reachable with the current X-ray and SZE observations, 
and equivalent for clusters of different redshift. One candidate is the virial radius. In a 
Friedman-Robertson- Walker universe, an unperturbed spherical region expands indefinitely, 
while a perturbed overdense region (the seed of a future cluster) eventually recollapses. 
When the overdense region collapses under the effect of its own gravity, it is ass umed to 
reach virial equilibrium when the radius is half of that at maximum expansion (IPeebles 
19801 ; lLacey fc Coldll993l ). The ratio of the mean cluster density to the background density 
at the time of virialization is A„ = 187r 2 for a universe with critic al matter density (Qm — !)• 



For a different cosmology with Qf, — 0, iBryan &: Norman! (119981 ) showed that A^ ~ 18 



82x - 39x 2 , where x = Q M0 {1 + z) 3 /E 2 (z) and E 2 (z) = fWl + z) 
found from a fit to numerical simulations (lLacey fc Coldll993l ). 



as 



With this characterization of the mean cluster density at time of virialization, the virial 
radius can be determined as the radius within which the average density of the cluster is A v 
times the critical density, via 
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irp c (z)A v (z)rl ir = M tot (r Av(z) ) 



in which both p c (z) and A v (z) are cosmology dependent, and the critical density p c (z) is 
defined as: 

3 V 



Pc(z) 



ZH 2 n E(z) 2 



8ttG 



(2) 



Unfortunately, the virial radius is usually unreachable with current X-ray and SZE 
measurements, and one is forced to perform measurements out to a smaller radius. Such 
a radius (r&) is characterized by the density contrast parameter A in place of A v (z) in 
Equation [1], and corresponds to a higher average density, 4/3 • irp c (z)A ■ r A = M tat (r A ). We 
choose a contrast parameter A = 2500, corresponding to an average density of 2500 times 
the critical density at the cluster's redshift. This choice is motivated by the fact that this 
is the radius typically reachable with our SZE and X-ray data without any extrapolation of 
the models (B2006, L2006).Q 



2.2. Scaling relations 



The hierarchical structure formation theory developed by iKaiserl (1l986l ) predicts simple 
relationships between physical parameters of collapsed structures, known as scaling relations. 
With the assumptions of hydrostatic eq uilibrium and of an isoth ermal distribution for both 
the dark matter and the cluster gas (e.g.. lBryan fc Normanlll998l ). it can be shown that there 



is a simple relationship between a cluster's total mass and its gas temperature T e : 



T P oc M 2 t l?E(z?/* 



L gas 



(3) 

tot — M tot (r A ). For 

M gas (r&)), the expected relationship between the gas mass 



where the mass is calculated out to a radius of mean overdensity A, M t 

fgas = M gas /M toU (M, 



within ta and the gas temperature is 



n2/3 



(4) 



1 The use of a constant overdensity factor A was shown by Maughan et al. ( 2006^ to give results similar 
to the case of a variable overdensity factor A(z) = A(0)[A K (^)/A tI (0)], in which the variable overdensity 
scales with redshift in order to keep the ratio of two comoving densities constant. 
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The Compton-y parameter is a measure of the pressure integrated along the line of sight: 



DC 



k B T e 

V= I cr T n e \dl (5) 



m e c 2 

One can further integrate the y parameter over the solid angle Q subtended by the cluster, 
to obtain the integrated Compton-y parameter: 

k B a T s 
' \ \m e c 2 j ,/„ ./.i 



Y = I y<M = W ( ) I dl j n e T e dA (6) 



where A is the area of the cluster in the plane of the sky. In the context of an isothermal 
model, Y is proportional to the integral of the electron density n e over a cylindrical volume, 
thus 

r B5 o,T,/n, ff = Vl = W . (7) 

In section we consider the effect of integrating gas mass within a spherical volume while 
determining Y in a cylinder. Using equation [3] we can rewrite equation [7] in terms of either 
M tot or T e , or substitute M gas / f gas for M tot , to obtain: 

YD\ oc f gas T^E{z)- 1 
YD\ oc f gas M?i 3 E(z) 2 / 3 (8) 

Equations M are the scaling relations that we investigate observationally in this paper. 



3. SZE and Chandra X-ray observations of galaxy clusters 

3.1. Data 

We analyze the SZE and X-ray data observations of 38 clusters in the redshift range 
z=0.14 — 0.89, observed with the BIMA and OVRO interferometric arrays and with the 
Chandra X-ray imaging spectrometers. Both data modeling with the isothermal f3 model 
and the data themselves are presented in B2006 and L2006, the previous two papers in 
this series. We re fer to L2006 for details on the observations and data modeling, and to 



Reese et al.l (120021 ) for a detailed illustration of the modeling of the OVRO/BIMA SZE data 
in the Fourier plane. In the following, we review those aspects of the data modeling and 
analysis that are relevant to the investigation of the scaling relations. 
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3.2. Data modeling 



The ga s density model is based on the spherical /3-model (ICavaliere fc Fusco-Femiano 



19761 . Il978l ) , which has the form 



n, 



1 + 



-3/3/2 



where n e o is the central electron number density, r is the radius from the center of the cluster, 
r c is a core radius, and (3 is a power-law index. When integrated along the line of sight to 
determine the projected SZE decrement distribution (oc n e ) and X-ray surface brightness 
(oc rig), this model has the simple analytic forms 



AT 



S , 



AT n 



02 

e 2 



(l-3/3)/2 



(l-6/3)/2 



(9) 



(10) 



where AT is the central thermodynamic SZE temperature decrement /increment an d 9 C is 



the angular core radius of the cluster (e.g., iBirkinshaw et al.lll99ll ; iReese et al.ll2002l ). This 



model typically provides a good description of the X-ray surface brightness and SZE decre- 
ment profiles out to ^r^nn fe.g.. IJones fc Formanlll984j ; lElbaz et al.lll995c iGrego et al.ll200ll ; 
Reese et al.lbood : lEttori et alJboO j i7 This simple model, however, does not provide a good 



description of the peaked X-ray surface brightness observed in the center of some clusters. 
To minimize the systematic bias associated with modeling of cluster cores, we therefore ex- 
clude the central 100 kpc from both the spatial and spectral X-ray data, as was done in the 
previous two papers in this series. The emission-weighted X-ray spectroscopic temperature 
is also determined by a single-temperature fit to the X-ray spectrum of photons extracted 
from an annulus between 100 kpc and r25oo (L2006). This 100 kpc-cut model was shown to 
recover the gas masses of simulated clusters with a range of dynamical states to better than 
5% accuracy at r 2 5oo (L2006) and, when applied to the determination of the Hubble constant, 
yielded the same results as a more complex non-isothermal model (B2006). Uncertainties 
associated with the isothermal assumption are included via an additional systematic error 
described in section H~T1 

Following our earlier analysis methods, we do not model the dark matter distribution, 
and calculate the total mass directly from the equation of hydrostatic equilibrium (see Equa- 
tion (TTJ below). An upcoming paper (Mroczkowski et al. 2008, in preparat ion) will extend 



the a nalysis to use more sensitive data from the Sunyaev-Zeldovich Array (IMuchovej et al. 



2007) and more a c curat e mo deling of the cluste r gas, based on the non-isothermal models 
of IVikhlinin et al.l (bood ) and iNaeai etaD J2007h . 
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3.3. Analysis Methods 

Best-fit model parameters and confidence intervals for all model param eters are obtained 



using a Markov chain Monte Carlo (MCMC) method described in detail by lBonamente et al. 



(120041 ) and B2006. L2006 explains the implementation of the likelihood calculation for the 
100 kpc-cut model. For each cluster, the Markov chain constrains the parameters S x q, ft, 9 C , 
AT , T e , and abundance (see L2006 for best-fit values). We use the cosmological parameters 
h = 0.7, Qm=0-3 a nd ^a=0.7 to calculate each cluster's angular diameter distance Da (e.g., 



Carroll et alJll992h . 



From these model par ameters we calcula te r 2 5oo and M tot {r 25 oo) through the hydrostatic 



equilibrium equation (e.g.. lGrego et al.ll200ll ). 



G/jrrip r 2 + r 2 ' 

in which /i is the mean molecular weight calculated using the X-ray metallicities and r c = 
c Da- One obtains r 2 5oo from the solution of the following equation: 

in which p c (z) is given by Equation [2J A = 2500, and the right hand side is just M iot (r 250 o). 
We then compute the global cluster quantities needed for the analysis of SZE and X-ray 
scaling relations. The gas mass is computed by integrating the gas density model, 

f 1 + 9 2 d6, (12) 

where fi e is the mean molecular weight of the electrons, and n e o is the central electron density, 
obtained from the parameters of the f3 model (L2006, Equation 12). 

The integrated y parameter (Y, equation E]) is calculated using the measured SZE decre- 
ment AT, which is directly proportional to the Compton-|/ parameter 

AT = Tomb • f(x) ■ y. 
The factor f(x) is the frequency dependence of the SZE: 

f(x) = (x ■ coth (|) - 4) ■ (1 + 5 rel ), 

in which x = hu/kBTcMB, and 5 re i is a small relativistic correction factor. At our observing 
frequencies, f(x) ~ —2. Thus, 

f AT AT rr2500 ID A / q2\ (l-3/3)/2 

y= / ^773^ = ^73 / 1 + 7^ on. (is) 



- o 

< MBf{xy U " TcMBf(x) Jo V ' °c 
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4. Observational constraints on SZE scaling relations 
4.1. Regression method 



Our measurements of masses and integrated y parameters, using the method described 
in Section [31 are shown in Table HJ The errors in Table [1] represent the photon-counting 
statistical uncertainties of the X-ray data and the statistical uncertainties of the SZE obser- 
vations. Additional sources of uncertainty in the measurement of cluster parameters include 
cluster asphericity and projection effects, small-scale clumping of the gas, the presence of 
point sources in the field, CMB an isotropy, the assum ption of isothermality, and instru- 
mental calibration, as discussed by iReese et al.l (120021 ). L2006, and B2006. Therefore, in 



fitting Y versus T x , M gas , and M tot (Equations [8]), we include an additional statistical error 
(combined in quadrature) of ±20% for the masses, and ±10% for Y and T x . 

We perform a line ar least-squares re gressi on in log spa c e, loq (Y) = A ± B ■ log(X), 
following the method of iPress et al.l (119921 ) and lBenson et al.l (120041 ) . This method accounts 
for errors in both measured parameters for each scaling relation, and it minimizes the \ 2 
statistic defined as 



X 



E 



(logjYj) -A- BlogQQf 
a l 9 (Y t ) + (Ba log(Xl )) 2 



in which cr log(Yx ) = a Yi /Y i log(e), v log{Xi) = a x J Xilog{e), and the linear errors a Yl and o Xi 
are obtained from the upper and lower uncertainties around the best-fit values as o = 
(a+ + a-)/2. 



4.2. The Y-M gas , Y-M tot and Y-kT scaling relations 



The above derivation of the self-similar scaling relations does not include any variation of 
the gas fraction with cluster m ass. However, there m ay be some eviden ce for such variation i n 
both X-ray observations (e.g.. IVikhlinin et al.ll2006l ) and simulations (jKravtsov et al.ll2005l ). 
We examine this in the present data by performing a logarithmic fit to the f gas -M gas data 
using a linear relationship (log(Y) — A + Blog(X)). We find no significant evidence for a 
variation of f gas with mass (B = 0.14 ±0.08). In the following we therefore assume that f gas 
is a constant. We then perform similar logarithmic fits to the Y-M gas , Y-M to t and Y-kT 
data. The results are shown in Figure [2] and in Table [2J Under the assumption of a constant 
fg as , all scaling relations are consistent (within 2a statistical uncertainty) with the simple 
self-similar model of cluster evolution. 
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4.3. Redshift evolution of the Y-M gas , Y-M tot and Y-kT scaling relations 

The large number of clusters (38) and redshift coverage of our Chandra and OVRO- 
BIMA data (0.14 < z < 0.89) enables the investigation of a possible redshift evolution of 
the SZE scaling relations. For this purpose, we divide our sample evenly into low-redshift 
clusters (z < 0.30, 19 clusters) and high-redshift clusters (0.30 < z < 0.89, 19 clusters), and 
repeat the logarithmic fits of Section 14.21 

The results of Table [2] indicate no evidence for redshift evolution of the SZE scaling 
relations, as f gas is consistent with a constant for both low and high-redshift clusters. Fur- 
thermore, the SZE scaling relations are consistent with the self-similar slopes at or below 
the 2.5cr level. 



5. Comparison of SZE and X-ray Measurements 

In the previous section we have examined the SZE scaling relations based on quantities 
derived jointly from SZE and X-ray observations. Here we compare the cluster properties 
derived from our SZE observations without using the X-ray data in the fits, in order to 
determine whether the relations we observe depend strongly on the X-ray information. 



We analyze the SZE data using the model of Section 13. 2\ using additional assumptions 
to provide the constraints that would otherwise be provided by the X-ray data. As a first 
assumption we fix (3 = 0.7 and fit for 8 C and ATq in Equation [9j following L2006. The choice 
of fixing (3 = 0.7 is determined by the fact that this is the median value for our sample; L2006 
also show that using values of 0.6 and 0.8 results in changes to the parameters that are small 
relative to the 68% statistical uncertainties. The data quality allows us to perform this SZE- 
only analysis for 25 of the clusters in the full sample, as shown in Table [3] Knowledge of the 
gas temperature is required in order to determine r25oo , as can be seen from the combination 
of Equations [1] and [TT] 



' ( 3(3kT \ 1 2 

' c ■ 



r250 ° V \GpmJ ln Pc (z) • 2500 



In the absence of complementary X-ray spectroscopic data, we estimate th e gas temperatur e 



directly from the SZE data following the iterative method described by I Joy et al.l (120011 ). 
We first choose an initial estimate of the gas temperature, from which we obtain r25oo and 
M tot (r 2 5oo)- We derive M 9as (r 2 5oo) using the equations described in Section [373| with the 
central gas density n e o calculated from the parameters of the SZE decrement model (L2006, 
Equation 13). We provide a final constraint by assuming that the gas mass fraction of each 
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cluster is equal to the average value for this sample, f gas = 0.116 (L2006), and iteratively 
lation M gas (r 2500 ) = f gas ■ M tot (r 2500 ) in ord 

T and r 25 oo- 



solve the equation M„ as {r 2 ^o) = f qas • M tot (r25oo) m order to find self-consistent estimates of 



Results of the SZE analysis are shown in Table [31 In Figure [3] we compare the SZE 
measurements of gas temperature, r25oo, Y, and gas mass with the values from the joint 
analysis of Section [3j In the joint analysis, the X-ray data are solely responsible for the 
measurement of the temperature, and drive the fit of the spatial parameters and gas density. 
They therefore drive the measurements of r 2 5oo and M gas as well. The quantities inferred 
from the SZE data are in good agreement with those from the joint analysis, indicating 
that we have not altered the scaling relations of Section 14.21 by incorporating the X-ray 
data. Although these results show that it is possible to estimate Y and the cluster gas mass 
from the SZE data alone, the joint analysis is preferred because it provides the strongest 
constraints on cluster properties, requires fewer assumptions about the cluster structure and 
composition, and can be applied to the full 38-cluster sample. Nevertheless, the ability to 
derive cluster properties fro m SZE observations will be important for the many SZE cluster 



surveys currently underway (iRuhl et al.ll2004J : lFowlerll2004j ; lKanekoll2006l ): these surveys are 
expected to generate large samples of SZE-selected clusters but will not generally have access 
to deep X-ray observations for cluster character ization. The gas te mperature may also be 



inferred from multi- frequency SZE observations ([Hansen et al.ll2002l ). 



Recent work by lKravtsov et all (120061 ) indicate that the quantity Y x = M gas ■ k B T e , the 



X-ray analogue of Y, is a low-scatter proxy for the cluster total mass. For the isothermal (3- 
model used in this paper, the gas mass can be estimated from X-ray data by using Equation 
T2l and thus our data also provide a measurement of Yx- We can therefore compare Y with 
Yx, in order to establish observationally whether the two quantities are indeed equivalent. 
In the case of the isothermal /3-model, the integrated Compton-?/ parameter is an integral of 
the electron density over a cylinder C of infinite length along the line of sight, and of area 



„2 . 
2500- 



n e (r)dV 



k B a T T e 



1 



n e m P ii e dV. 



m e c* J m p /i e j c 

Since the gas mass is given by an integral over a sphere S of radius r 2 5oo, 



.1/ 



n,, m pi i, (IV 



the relationship between Y and Yx is 



YD\ 



m e c 2 



m p /j, e 



-CY X 



(14) 
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where the constant C = j c n e dV/ J s n e dV accounts for the different domain of integration 
of Y and Y x , and depends on the parameters of the (3 model. We calculate C separately for 
each cluster, and typically find C ~2. 

In Figure H] we plot Y, as derived from the joint SZE/X-ray analysis, against Yx- A fit 
of the data to the relationship in Equation [14] (the dotted line in Figure HI with no degrees of 
freedom) results in an acceptable \ 2 statistic, corresponding to a null hypothesis probability 
of 80.4%. The agreement with the relationship in Equation [T41 shows that Yx is an unbiased 
estimator of the integrated Compton-?/ parameter, within the uncertainties of the current 
measurements. 



6. Comparison with theoretical simulations 



We compare our re sults with those of recent cosmological cluster simulations (INagai 



20061 ; iNagai et al.l 120071 ) that include radiative cooling, UV heating, star formation, and 
stellar feedback processes in addition to the standard gas dynamics. In Figure [5] we compare 
16 clusters simulated at z=0 and 0.6 using cooling and star formation feedback processes (in 
red), the same sets of clusters performed using non-radiative gas dynamics (in green), and 
our 38 clusters observed with Chandra and OVRO-BIMA (in black). 

The best-fit power-law models that describe the two simulations are shown as dashed 
lines in Figure [51 A fit of our data to the cooling and star formation model (red dashed 
line) results in a x 2 null hypothesis probability of 99.9%, and the non-radiative model (green 
dashed line) has a probability of 0.5%. The comparison indicates that both simulation 
models show a similar slope to the observed clusters, with the cooling and star formation 
feedback model providing a better match to the data. 



7. Discussion and conclusions 



We have investigated scaling relations between the integrated Compton-?/ parameter Y 
and total mass, gas mass, and gas temperature using 38 clusters observed with Chandra and 
OVRO-BIMA. Fits of the Y-M gas , Y-M tot and Y-kT data to a power-law model agree with 
the slope predicted by a self-similar model in which the evolution of clusters is dominated by 
gravitational processes. The normalization o f the Y-M gas scaling relation agrees well with the 
numerical simulations of INagai et al.l (120071 ). in which collisionless dynamics of dark matter 
and gas dynamics are complemented by cooling, star formation and feedback phenomena. 
The agreement provides observational evidence that non-gravitational phenomena may also 
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be an important factor in the physics of clusters. 

The redshift coverage of our sample enabled an analysis of the scaling relations as 
function of redshift, by defining a low redshift sample (0.14 < z < 0.30, 19 clusters) and a 
high-redshift sample (0.30 < z < 0.89, 19 clusters). Both samples follow the prediction based 
on the self-similar model. Our data indicate no significant evolution in the SZE properties 
of clusters at redshift z < 1. 

We also measure the cluster mass and integrated Y parameter using the SZE data alone, 
without making use of the Chandra spectral and spatial information. These measurements 
are in good agreement with those based on the joint X-ray/SZE analysis, providing evidence 
that SZE surveys can be used to determine the number density of galaxy clusters as functions 
of mass and redshift. 
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Table I. Cluster Parameters from Joint Analysis of X-ray and SZE Data 



Cluster 


z 


D A E(z) 


»"2500 


kT 


M ga3 Mtot 


Y 






(Gpc) 


(") (kpc) 


(kcV) 


(10 13 M Q ) (1O 14 M ) 


(mjy) 



0.14 


0.52 


1.07 


0.18 


0.63 


1.09 


0.25 


0.81 


1.13 


0.17 


0.60 


1.09 


0.32 


0.96 


1.18 


0.23 


0.76 


1.12 


0.20 


0.68 


1.10 


0.15 


0.54 


1.08 


0.18 


0.63 


1.09 


0.16 


0.57 


1.08 


0.22 


0.73 


1.12 


0.23 


0.76 


1.12 


0.38 


1.07 


1.22 


0.17 


0.60 


1.09 


0.29 


0.90 


1.16 


0.18 


0.63 


1.09 


0.26 


0.83 


1.14 


0.28 


0.88 


1.15 


0.22 


0.73 


1.12 


0.54 


1.31 


1.34 


0.89 


1.60 


1.65 


0.58 


1.36 


1.37 


0.69 


1.47 


1.47 


0.54 


1.31 


1.34 


0.49 


1.25 


1.30 


0.55 


1.32 


1.35 


0.57 


1.35 


1.36 


0.48 


1.23 


1.29 


0.41 


1.12 


1.24 


0.55 


1.32 


1.35 


0.83 


1.57 


1.59 


0.78 


1.54 


1.55 


0.33 


0.98 


1.19 


0.58 


1.36 


1.37 


0.45 


1.19 


1.27 


0.81 


1.56 


1.57 


0.24 


0.78 


1.13 


0.29 


0.90 


1.16 



Abcll 1413 

Abell 1689 

Abell 1835 

Abell 1914 

Abell 1995 

Abell 2111 

Abell 2163 

Abell 2204 

Abell 2218 

Abell 2259 

Abell 2261 

Abell 267 

Abell 370 

Abell 586 

Abell 611 

Abell 665 

Abell 68 

Abell 697 

Abell 773 

CL J0016+1609 

CL J1226+3332 

MACS J0647.7+7015 
MACS J0744.8+3927 
MACS J1149.5+2223 
MACS J1311. 0-0310. 
MACS J1423.8+2404 
MACS J2129.4-0741 
MACS J2214.9-1359 
MACS J2228.5+2036 
MS 0451.6-0305 . 
MS 1054.5-0321 . 
MS 1137.5+6625 
MS 1358.4+6245 
MS 2053.7-0449 . 
RX J1347.5-1145 
RX J1716.4+6708 
RX J2129.7+0005 
ZW 3146 



e 



206+5 
219+P 
172±| 
228+^ 
133±| 
141+* 
206±| 
256+];? 
191±| 
172±| 
148+ 
132+ 
97±| 
182+1 
111±| 
162+f 
153+1° 
134+5 
148+f 
80±| 
66±| 
92±| 
59±| 
71±| 
74±| 
66+ 
73+2 

9i±- r ; 

81±| 
82±| 
89+1 

42 ±3 

113+J 

54± 5 
122±| 

45+4 
128+1 
132+3 



664+lf 

672+20 

660+JI 

621+21 

518+S 

682+10 

671+1 

581±i| 

476+S 

525+11 

484+11 

508±|J 

529±|g 

482±i° 

490+1 

616+1? 

568+2] 

527+fO 

507±i| 

512+ 

606+15 

420+g 

449+1 

448±|« 

422+H 

474+30 

547+22 

444+ 

526+JI 

686±g| 

311+H 

539+1? 

358+fl 

706±| 

341+1 

486+J2 

574+H 



8:? 



7.5±!!;t 
10.5+" | 

9.5±°;f 
8 2+0' 4 

8 - 2± rt 8 
14.8±0;| 

7.8±0;| 
5.8±»;| 
7.4+ 
5.9±, 
8-7±S| 
6.4±°;| 
6.8±»;J 

9.6+1;* 
in 2+° ^ 

8-2±°,| 
10.5±o| 
13-5+H 

i4.i± : 

8.1±«;| 
9.9±»;f 

7 - 2± d 

7 n+0'4 

8.6± 08 
10.2±i,-0 
8.4+81 



9.9±0;» 
9.8+F 
4.5+ 
8.9+ 

4 -8±o Y 6 
16.5+*;° 

6.6±i;i 

8.3+n'-) 



8:§ 
8:1 
8:? 



2.6±o 

5.1+ 

5.8+0 

4.8+0 

3.5+0 

2.2+0 

8.1+ 

4.7+0 

3.o±; 

1.8+ 
3.0+ 

2.2±: 

2.8+ 

2.3±: 
2.4±: 

2.6±; 
3.6+ 

4.4±: 
2.7±: 

4.4+ 
3.9+ 
4.9+' 
3.1+: 

3.i±: 

2.1+ 
2.3+ 
3.3+ 
3.9+ 
2.8±: 
4.8+ 
7.4+ 
1.2+ 
2.5+ 
0.9+ 
8.8±; 
1.2+ 

2.6±; 
4.4±: 



8:1 



i!:± 



[]■}, 



o:¥ 
1 
:i 
\ 



8:1 



2 2+°' 2 
5.0+0-f 
5.6+ 
4.8+ 
4.7+ 
2.5+ 
5.5+ 
5.0+ 
3.3+ 
1.8+ 
2.6+ 
2.0+ 
2.8+ 
2.5+ 
2.1+ 
2.0+ 
4.3+ 
3.5+ 
2.6+ 
3.3+ 
5.2+ 
6.0+ 
2.3+ 
2.3+ 
2.2+ 
1.9+ 
2.8+ 
3.9+ 
2.0+ 
3.8+^ 
4.5+ 
l.Oijj 
3.1+ 
1.2+ 
8.1+ 
1.4+ 
2.1+ 
3.6+0 



1:1 
4:1 



2 
9:1 
8 

2 



2.99+0-44 
3.79+pl 
2.09±«;i| 

0.95+ ;| 
6.89±°- BS 



§8 

m 



4.43+ 
1.94+ 
0.82+ 
1.34+ 
0.72+ 
0.71+ 
1.03+ 
0.54+ 
2.68+ 
1.01+ 
1.67+ 
1.68+ 
0.73+ 
0.35+ 
0.62+ 
0.34+ 
0.58+ 
0.28+ 
0.28+ 
0.39+ 
0.77+ 
0.94+ 
0.66+ 

0-77±n 1( , 
n n9+o-()i 

0.56+0 08 
0.09±°;| 
1.62+0-ls 



8:?I 

m 

81 

8:?e 

IS 
8:18 
8:A8 
8:88 



M 
:8§ 
■flfi 



0.10+ 
0.66+ 
0.88+ 



8:°? 
8:15 



0.120+8 
0.102+0 
0.104+0 
0.101±[j 
0.074+0 
0.088+2 
0.147+2 
0.096±g 
0.090+; 
0.101+ 
0.119+ 

0.H0+: 

0.100+ 
0.091+ 

0.1H+: 
o.i3i±: 

0.084+ 
0.126+: 

o.io6±: 

0.131+ 
0.075+ 
0.082+: 
0.136+! 
0.134+ 1 
0.097+ 

0.H6+: 

0.116+ 
0.102+ 
0.138+: 
0.128+ 
0.164+ 
0.115+ 
0.081+ 
0.076+ 
0.109+: 
0.088+: 
0.124+; 
0.122+ 1 



8:88? 
8:881 
8:89! 
8:8ftJ 



8:898 



8:8Ai 
8:898 



8:898 



Table 2. Scaling Relations from Joint Analysis of X-ray and SZE Data 



Scaling 


All clusters 


0.14< z 


<0.30 


0.3K z 


<0.89 


relation 














A 


B 


A 


B 


A 


B 


f gas gas 


-2.86±1.09 


0.14±0.08 


-2.60±1.79 


0.12±0.13 


-3.00±1.37 


0.15±0.10 


Y,kT 


-6.24±0.22 


2.37±0.23 


-6.33±0.32 


2.46±0.34 


-6.13±0.30 


2.27±0.30 


Y,Mgas 


-23.25±1.77 


1.41±0.13 


-25.86±3.45 


1.60±0.25 


-21.43±3.00 


1.28±0.15 


Y,M tot 


-28.23±3.00 


1.66±0.20 


-31.20±5.35 


1.87±0.35 


-25.45±3.46 


1.47±0.23 
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Table 3. Cluster Parameters from Analysis of SZE Data 



Cluster 


z 


' ZOUU 


kT 


M„„e 

gas 


- LV - L tot 


Y 






V ) 


(keV) 


(10 13 M n ) 


(10 14 M n ) 




Abell 1689 


0.18 


196±| 


8.0+°;? 


4.2±g : g 


3.6+8;! 


1-88+8:11 


Abell 1835 


0.25 


169±| 


"■4±5:§ 


fi 1 4-0.6 


5-3±8:I 


2.66+°;°° 


Abell 1914 


0.17 


212+ 1 , 


8.5±i;° 


4 c ,+°- 7 


3.9±8| 


2.39+8:11 


Al 11 1 Afl r" 

Abell 1995 


0.32 


117+3 




3-7±8:i 


3-2±o^3 


c\ o a i 16 


Abell 2111 


0.23 


124±| 


5-4±i;g 


2.0±8;I 


i.7±8:f 


0.45±8;| 


Abell 2163 


0.20 


229+1?, 


15.6+2-4 


8-7±i;l 


7.5±{i 


8.03± 2 ; 26 


Abell 2218 


0.18 


206+}? 


10.0+11 


4.8+°;° 


4 1+°- 8 


Q on_i_1.24 


Abell 2261 


0.22 


146±| 


6.4±i;g 


2-9±8;g 


z - o=t 0.4 


0-76+8:11 


Abell 267 


0.23 


138±£ 


6.5±i;° 


2-7+8:1 


2.3±8;| 


o.75±8:i? 


Abell 370 


0.38 


96±| 


7.3±8;f 


3.o±8:i 


2-6±8:i 


o.45±8:iB 


Abell 665 


0.18 


181±i? 


1A±\\ 


o Q+0.7 


2-8±8| 


1.55+8;!? 


Abell 697 


0.28 


130+J3 


8.3±?| 


o 7-1-1.2 
0.8 


3.2± ; 7 


o.99±°;i 


Abell 773 


0.22 


150+°, 


6.9±i;° 


3.i±8:i 


2.7±8| 


o.93±8:I 


CL J0016+1609 


0.54 


81±| 


12.6+°;° 


4-i±8:t 


3-5±o:4 


0.92+8;}} 


CL J1226+3332 


0.89 


53+2" 


9-2±8:l 


3 2+°- 3 


2-7±8:i 


o.27±8;| 


MACS J0647.7+7015 


0.58 


72±| 


8.6±i;i 


3.4±°;| 


2-9±8:i 


0.38± ; 09 


MACS J131 1.0-0310. 


0.49 


73±| 


6.8±°;I 


2-4±8i 


2.i±8:l 


o.26±8:8? 


MACS J2214.9-1359. 


0.48 


91±i 


10.2±°-| 


4.5±°;| 


3.9±8:i 


o.74±8;l| 


MACS J2228. 5+2036 


0.41 


100+2 


li.o+ii 


4-3±8;| 


o 7_|_0.4 
■ 0.4 


i.os±8:i? 


MS 0451.6-0305 


0.55 


84±| 


11.3±Ji§ 


A 7+ 0A 


4.i±8:l 


0.79±8;?| 


MS 1137.5+6625.... 


0.78 


49±| 


6.4±°i 


i-9±8:§ 


i-6±8:§ 


o.i2±8;8l 


MS 1358.4+6245.... 


0.33 


99±| 


6.4+°;° 


2-4±8;| 


2.i±8:i 


0.40±8;ol 


RX J1347.5-1145.... 


0.45 


108±| 


i2.7±l:5 


6.5±°; 6 7 


5.6+°;° 


1.41±8;i 


RX J2129. 7+0005... 


0.24 


123±! 4 


7.9+1;? 


2.i±S:B 


i-8±8:l 


o.80±8;i 


ZW 3146 


0.29 


128±g 


8.3± ; 9 


3.9±°;| 


3.3+8:1 


o.94±8:i 



o 



HEIH 



PO 

X 

LO 



o 

X 



© 



O 



Z'O 



L'O 



sd6 



ceo 



Fig. 1. — Dependence of / ffas on M 9as ; open squares (in black) are clusters at 0.14 < z < 0.30, 
open diamonds (in red) are clusters at 0.30 < z < 0.89. The gas fraction at r 2 5oo shows no 
evidence of evolution with mass for the clusters in this sample. 
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Fig. 2. — Scaling relations between the SZE Y and M gas , M tot and kT e . Open squares (in 
black) are clusters at 0.14 < z < 0.30, open diamonds (in red) are clusters at 0.30 < z < 0.89. 
All measurements follow simple power-law models with indices that are consistent with the 
values of the self-similar scaling theory (Tabled]). 
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Fig. 3. — Comparison between cluster parameters derived from joint analysis of X-ray and 
SZE data (x axis) and those derived from SZE data (y axis) following the procedure of 
Section Dashed lines correspond to y = x. 
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Fig. 4. — Comparison between the SZE Y parameter and the X-ray quantity Yx (see Section 
[5] for explanation of the normalization constant); open squares (in black) are clusters at 
0.14 < z < 0.30, open diamonds (in red) are clusters at 0.30 < z < 0.89. The dashed black 
line corresponds to y = x. 
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Fig. 5. — M 9as vs. F for simulated and observed clusters. Open squares (in black) are 
OVRO/BIM A/Chandra measurements. Also shown are simulated clusters from a cooling 
and star-formation feedback model (red c ircles) and simulated clusters from a non-radiative 
model (green triangles) from iNagail (120061 ). For comparison, all of the M gas and Y quantities 
are integrated over a spherical volume, as described in Section [5j Open symbols represent 
simulated clusters at z=0 and filled symbols represent simulated clusters at ,2=0.6. The 
power-law fits to the simulated clusters (dashed lines) are A = — 25.83±0.52, B = 1.58±0.04 
for the cooling/star-formation model and A = —25.79 ± 0.80, B = 1.56 ± 0.06 for the non- 
radiative model. 



